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We show that an extra constant of motion with an analytic form can exist in the neighborhood of 
some discrete Wannier orbits of helium, when one includes retardation and self-interaction effects. 
The energies of these discrete stable Wannier orbits are in the correct atomic magnitude. The 
highest frequency in the stable manifold of one such orbit agrees with the highest frequency sharp 
line of parahelium to within 2%. The generic term of the frequency in the stable manifold to higher 
orbits is also in agreement with the asymptotic form of quantum mechanics for helium. 

PACS numbers: 31.15.Ct, 03.20.+i,05.45.+b 



In this letter we explore some surprising consequences 
of the retardation effects of Maxwell's electrodynamics 
to atomic physics. We show that electrodynamics with 
retardation prescribes a discrete set of stable Wannier 
orbits for the helium atom. A Wannier orbit is one where 
the two electrons are in the same circular orbit and in 
phase opposition, that is, along a diameter Jll. 

The linearized dynamics about a Wannier orbit has 
a one-dimensional stable direction, a one-dimensional 
unstable direction and ten neutrally stable oscillatory 
modes For an infinitely massive cv-particle, helium 
is a six-degree of freedom Hamiltonian system with only 
four independent constants of motion, (namely, the en- 
ergy and the three components of the total angular mo- 
mentum). In the neighborhood of Wannier orbits, even 
if one restricts to the center manifold one still has 
a five-degree of freedom system. Because there are only 
four constants, the dynamics in the neighborhood of a 
generic Wannier orbit can diffuse away to large distances 
of the orbit. Here we show that an extra complex con- 
stant can exist in the neighborhood of a discrete set of 
Wannier orbits, when one includes retardation effects. 

To include the retardation, we use the work of Page 
(1918) || on the expansion of the Lienard-Wiechert for- 
mula p|. This formula is complicated because of the 
retardation constraint, and one way to convert it into a 
useful differential equation is to develop the constraint in 
a Taylor series. This was done by Page up to the fifth 
order in 1918 ||. Truncated to second order in (v/c), 
this retarded interaction is described by the Darwin La- 
grangian (J], which introduces the first retardation cor- 
rection to the Coulomb interaction. 

The Darwin Lagrangian is rotationally invariant, 
which generates an angular momentum-like constant of 
the motion according to Noether's theorem J(|. This con- 
stant is the angular momentum of the electrons plus the 
small correction of the angular momentum in the elec- 
tromagnetic field. There is also an energy-like constant 
(because the Darwin Lagrangian is time independent). 
According to a recent result, the existence of a constant 
of motion analytic in the neighborhood of an equilibrium 
point requires the eigenvalues of the linearized dynam- 
ics to satisfy some resonance condition Q. Since the 
eigenvalues in the presence of retardation depend on the 
particular orbit, an extra complex constant is possible for 
some discrete orbits, which is a genuine signature of the 



nonlinear dynamics with retardation. We stress that it is 
absolutely essential to include the retardation to unfold a 
degeneracy of the Coulomb dynamics, and not to attain 
some better precision. Last, we consider only the dynam- 
ics of spinless pointlike charges, and issues like quantum 
mechanical spin are not discussed. 

In this work we assume that the self-interaction ef- 
fects are described by the Lorentz-Dirac equation (LDE) 
with a renormalized mass [f|j|] ■ The LDE equation for an 
electron of charge — e can be written in the convenient, 
non-covariant form as 
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where x e is the electron's velocity, M e the renormalized 
electronic mass [|| and 7=1/ ^/l — (|x e |/c) 2 . In (Q) , 
F ext is the external electromagnetic force and in this 
work we use F ra d = 3^3 x e , which is the lowest order 
approximation in powers of (v /c) to the exact relativistic 
radiaton force (8). For a circular orbit, this is a term of 
order (v/c) 3 . 

We now introduce the expansion of the retardation 
constraint of the Lienard-Wiechert interaction, hence- 
forth called the Page series Let xi be the position 
of electron 1, of charge — e, and (3\ its velocity vector di- 
vided by c. The formula for the electric field made by 
electron 1 at a position X2 is 
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and 
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(3), 
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where n is the unit vector in the direction of X2 — xi and 
R12 = |x2— Xi|. Notice that all functions are evaluated at 
present time. For a circular orbit, in units of the zeroth 
order Coulomb term, (||) is of order {v/c) 2 and (^) is of 
order {v/c) 3 . We assume that the Taylor series of (||) 
converges, since (v/c) <C 1 in atomic physics. The first 
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term of the series for the magnetic field of electron 1 at 
x 2 is i 



B^ ) (x 2 ,xi,xi) 



2 [Pi x n] 
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In this work we consider only the above terms of the elec- 
tromagnetic interaction, as well as the relativistic correc- 
tion to the electronic dynamics up to order (v/c) 2 . 

Along a Wannier orbit, the total force acting on elec- 
tron 2 can be calculated, using (||), (||), (Eh and the self- 
interaction of electron 2, to be 



g<-x I+ *,)-70(l->S), 
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where R is the radius of the Wannier orbit. In (||) we 
have also added the Coulomb attraction of the a-particle 
(of charge 2e), that we suppose infinitely massive and 
resting at the origin. Along the Wannier orbit, because 
the electrons are in phase opposition, the force along the 
velocity cancels out, demonstrating that the Wannier or- 
bit is a possible periodic solution of the electromagnetic 
equations up to third order in (v/c). 

Notice the appearance of the dipole term in (O) , 



D = -e(ki +x 2 ). 
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The total far field caused by the electrons depends lin- 
early on the quantity D of (Q), up to quadrupole terms 
j| . If D is zero, the orbit does not radiate in dipole. The 
fifth-order quadrupole terms would be important only 
in a much longer time scale, of order T/(v/c) 5 In 
this letter we consider only the dynamical system de- 
fined by truncating the Schott-Page interaction and self- 
interaction to third order in (v/c). 

Next we consider the dynamics of helium in the neigh- 
borhood of a Wannier orbit, according to the Darwin 
Lagrangian. For the Coulomb dynamics, this stability 
roblem has already been considered by many authors 
~(|. Here we consider only the second order of the 
Page series, which is essential to unfold a degeneracy of 
the Coulomb dynamics. The detailed consideration of 
the dissipative effect of the third-order term of the Page 
series and the (also third-order) self-interaction term will 
be published elsewhere ||. Here it suffices to say that 
these higher-order terms are important only in a much 
longer time-scale. To second order in (v/c), the equations 
of motion for helium with retardation and relativistic cor- 
rection are 

d 1 2e 2 e 2 

Me ^ [(1 + 2 |il|2)il] = ~ ~R? (xi) ~ Bg^ 2 ~ Xl) 

+E 2 2) (xi,x 2 ,x 2 ,x 2 ) - -xi x B 2 ^(xi,X2,x 2 ), 
c 

d 1 2e 2 e 2 

Me ^ [(1 + 2 |X2 ' 2)X2] = ~ w (X2) ~ w^ 1 ~ X2) 

+Ei 2 J (x 2 ,Xi,xi,xi) x 2 x Bi 2 i (x 2 ,xi,xi), (8) 



where M e is the electronic mass, R\ = |xi|, i? 2 = |x 2 |, 
R\2 = |xi — x 2 |, E^l (x 2 ,xi,xi,xi) is defined in (0) and 
B^ 2 ' ) (x 2 , xi, xi) is defined in (||). The Wannier orbit is 
the following periodic solution of (||) 

x\ = Rcos(ujt), yi = Rsm(ujt) , z\ = 0, 
X2 = — -Rcos(wi), y 2 = — Rsm(ujt), z 2 = 0. (9) 

According to (||), the frequency of the orbit is related to 
i?by 
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which is correct to second order in (v/c). 

Linearizing (^J) about the Wannier orbit ([)[), we obtain 
a parametric linear differential equation with coefficients 
periodic in time and period T — tt/lo ||. If the Floquet 
exponents are all nondegenerate, one can find a complete 
set of solutions of the form JTl| l 

<5xi = exp(2zcij/ii) x^ exp(2inujt), 

Sx.2 = exp(2zw/ii) x 2 exp(2inuit), (11) 

n 

where the Floquet exponent /i is a complex number de- 
fined in the first Brillouin zone, —1/2 < Re(fi) < 1/2. 
Notice that for the Floquet components we use an upper 
index, but to label coordinates as functions of time, as 
in (]|), we use a lower index (to distinguish it from the 
Floquet components). To bring the variational equations 
to normal form, we define the coordinates 
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Next we calculate Hill's secular determinant [0 for the 
linearization of (||) about (||). As a simplification, let 
us define n = n + \i to be the running variable in the 
summations of ([TT|). Last, to introduce physical intuition, 
it is convenient to define the radiation coordinate x r and 
relative coordinate x,i as 



x r = X! + x 2 , 
x d = xi - x 2 . 



(13) 
(14) 



It turns out, as one can easily verify, that the linearized 
dynamics for the 5xd variable is decoupled from the 5x r 
variation. Here, for lack of space, we only develop the dy- 
namics of the Sx-d variable along the z = plane, correct 
to second order in (v /c) . We refer to [jlj for the Coulomb 
eigenvalues of the <5x r dynamics. A complete study of 
the eigenvalues will be published elsewhere ||. Before 
we write the planar variational equations, we need still 
another definition 
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We make the quantities U£ and V^f into the x and y 
components of a bidimensional vector K£ = (U^,V^), 
and also define K£ = (— V£,U£), x n = (x n ,y n ) and 
x„ = (— y*,x%). The variational equations for planar 
motion along the variable 5x,j can be written most simply 
in Floquet components as M 
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(x^ + 3K^ 
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We recall that n = n + fi, and there is one equation 
for every value of n. In principle, this would lead to an 
infinite Hill's secular determinant [0. For the case of 
circular orbits, when we use the £ and x variables, we 
find that the variables couple only to Xn+i an d vice 
versa. The variational problem for the coordinate Jx^, 
as defined by (16), turns into the following 2x2 matrix 
equation 
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where P 2 (n) = (-Ml + « + Q 2 (n) = {^f- - ^ ■ 
t^) and n + = (1 + n). We now consider the secular 
equation associated with jT7|): At = 0, we obtain the 
eigenvalues of the Coulomb stability problem 



(18) 



n = 0,-1,-1/2,-1/2. 



The next step is to consider the secular equation of 
for finite values of \(5\. The essential ingredient for us 
here is the degenerate root n = —1/2 of (|l8|). Because 
of this degeneracy, the double root n = —1/2 acquires 
a correction proportional to \f3\, instead of the generic 
second-order correction of the other nondegenerate roots. 
Solving the secular equation of ([l7]) , in the neighborhood 
of h — —1/2 and in the lowest order approximation in \f3\, 
we find 



-l/2± 



(19) 



According to ( pT| ) , h is a frequency of oscillation in units 
of 2lo. For example, along the linear mode of (|l7|) defined 
by n — —1/2, the frequencies in ( |Tl| ) are —to — 2ton and 
uj = 2a; (1 + n). 

Furta's theorem gives the condition for the existence of 
an analytic constant in the neighborhood of a fixed point 
of an autonomous vector field fh . The Wannier orbit is a 



fixed point of the autonomous vector field describing the 
dynamics in the system rotating with the frequency of 
the Wannier orbit. To apply the theorem, we must move 
to this rotating system. In this rotating system, the new 
frequency for oscillation along the mode (|l^) is found 
by adding 1/2 to ([l9|). In units of 2w, the frequency of 
Coulombic oscillations along the z direction in the mode 
Sz r is found in [pi to be 



- ~ 0.5345 . 



(20) 



and for planar Coulombic oscillations along the mode 
<5x r , in the rotating frame, is jl| 



LJ 2 




0.5560. 



(21) 



also in units of 2uj. A new resonance condition involv- 
ing the root ( |l9| ) is the easiest to be satisfied for small 
values of \j3\. Therefore, we are looking for a resonance 
among the frequencies of (pl|), (20), and (|l9|), which in 
the rotating frame evaluates to 



1.34961 



(22) 



in units of 2lo. According to Furta's theorem Q|, for 
an additional analytic constant to exist one needs a new 
resonance among the linear frequencies. By inspection, 
we find that the resonance involving the above three fre- 
quencies and with the minimal integer multipliers is of 
type 



LO\ — Ll>2 + 2nLL>3 = 0, 



(23) 



where n is an integer number and the factor of two mul- 
tiplying LU3 is because mode 3 comes from a degenerate 
eigenvalue. Condition (E3I) is satisfied for given by 



0.007969423 



(24) 



The binding energies, in atomic units, for the above val- 
ues of 1/31 are 



E = 



2.3853 



n~ 



(25) 



and the frequencies of the special "stable" Wannier orbits 
described by (^3|), in atomic units, are given by 



0.7442 



(26) 



Because we are involving a resonance among the planar 
oscillations 5x r , Sxa, and the 6z r oscillation, the lead- 
ing term in the McLaurin series [j?J of the constant is a 
product of Sz r times the planar coordinate variations. In 
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reference || we show how to construct this extra constant 
and discuss the shape of its level surfaces. 

The frequency lu z for oscillations along the 8z r linear 
mode can be obtained, as prescribed by (pp|), by multi- 
plying the above frequency by 2W2/7, 



0.7956 



(27) 



The inclusion of the third-order radiative terms makes 
this linear mode of ( p0| ) a stable mode ||, whereby the 
perturbation radiates energy and decays to the perfectly 
circular orbit. The frequency of this stable oscillation 
should then appear in the spectrum of sharp lines, given 
that the dynamics selects the neighborhood of the stable 
orbits. For example, the transition 2 1 P — > l^S of para- 
helium corresponds to a frequency of 2.9037 — 2.1237 = 
0.7799 atomic units ||Jl2|, and @ with n = 1 eval- 
uates to 0.7956, which is a 2% difference. The transi- 
tion 3 1 / 7 — * 2 1 5 in parahelium Jjl| has a frequency of 
2.1459 - 2.0551 = 0.0908 atomic units, and (jfj]) with 
n = 2 evaluates to 0.0995, which is an 8% difference. The 
asymptotic form of the quantum mechanical energies of 
helium are pi 



(z-i) 2 



2(n r + L + l) 2 ' 



(28) 



with Z = 2. For transitions with AL — 1, the frequencies 
predicted from the above formula are given by 



1 



L + l) 



(29) 



which agrees with ( p7| ) to within 20%. We stress that 
we relied only on dynamics to calculate the frequency of 
(p7|). Equation (|2^), on the other hand, was obtained 
in the usual way of taking differences between quantum 
energies, and was included only for comparison. 

A more detailed comparison with the spectrum of he- 
lium requires the consideration of elliptical orbits and 
zero-angular- momentum orbits ||. In reference 0] we 
include the third-order terms and find that the resulting 
damping of the oscillation mode of (|27j) agrees well with 
the quantum calculations of the linewidth. 

The extra complex constant is actually two real func- 
tions of phase space . Together with the four global 
constants of helium, there are six functions that can pro- 
duce a barrier in phase space, at least in the small region 
inside the resonance islands H]. Projected on the plane 
of the stable and unstable Wannier directions, the one- 
period return map of the full dynamics has a hyperbolic 
fixed point. The resonance islands are possibly found in- 
side the homoclinic tangle of the invariant manifolds, but 
this remains to be investigated in detail. 

Historically, the Wannier instability, first found by 
Nicholson pO] , led Bohr to postulate that, along a quan- 
tized set of orbits, some special mechanism would super- 
sede the Wannier instability, thus stabilizing the atom, 



which proved to be a very fruitful intuition Jl3| . Some 
of the details of how the six constants stabilize the dy- 
namics, in a time scale of order T/(v/c) 2 , are discussed 
elsewhere ||. 

The existence of this extra constant is a genuine ef- 
fect of electrodynamics with retardation. It is surprising 
that classical electrodynamics with retardation can offer 
such a sensible approximation to the spectrum of sharp 
lines. We do not know of any prior similar application of 
dynamics to atomic physics. 
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